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Abstract 

We analyse the abelian sandpile model on for the starting con- 
figuration of n particles in the origin and 2d — 2 particles otherwise. 
We give a new short proof of the theorem of Fey, Levine and Peres |3] 
that the radius of the toppled cluster of this configuration is 0{n^^'^). 

Consider the grid with an integer written inside each cell. We call this 
a configuration and think of the numbers as particles. We want to study a 
one player game with the following rule: whenever we see at least 4 particles 
in a cell {x,y), we are allowed to topple this cell, i.e. to remove 4 particles 
from {x,y) and add one particle to each of its neighbour cells {x,y — 1), 
{x,y + 1), (x — 1,1/), {x + 1,1/). The order in which we topple active cells, 
that is cells with at least 4 particles, is subject to our choice. Once no cell 
contains more than 3 particles, the procedure stops. Cells with less than 4 
particles are called stable and the whole procedure is called stabilisation. 

What we just described is called the abelian sandpile model. It is an exam- 
ple of a cellular automaton with a particularly simple 'symmetric' transition 
rule. The sandpile model was first defined in a book for school teachers in 
1975 by Engel 0, who called it the chip firing game. It was re-discovered 
in 1989 by the physicists Bak, Tang and Wiesenfeld [T] as a model for grav- 
itation (hence the current name). Later the sandpile model proved to be an 
extremely useful mathematical object, having connections, besides physics, 
to group theory, probability theory and partial differential equations. To get 
a feeling for the sandpile and to enjoy its beauty we refer the reader to the 
wonderful sandpile applet by Maslov The reader is welcome (and urged!) 
to check the proofs using Maslov' s applet. 

The sandpile model can be defined on a general graph by assigning par- 
ticles to the vertices and toppling a vertex whenever the number of particles 
assigned to it exceeds the degree of the vertex. In this text, however, we deal 
only with 'Z'^, d > 2. We shall use d = 2 for simplicity, however our results 
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generalise easily to d dimensions. We also, purely for aesthetic reasons, imag- 
ine the particles living in the cells rather than on the knots. Consequently, 
we distinguish one cell to be the origin of a cartesian coordinate system and 
assign coordinates to each cell. 

The word 'abelian' stands for the so-called abelian property of the sand- 
pile model: if the stabilisation terminates for some order of topplings, it 
terminates for any order of topplings. Moreover, for any two orders of top- 
plings that terminate the procedure, the toppled cells are the same up to a 
permutation and, consequently, the final configurations are identical. This 
allows us to speak of the final configuration resulting from some initial con- 
figuration. The proof of this statement is fairly standard and can be found 
in many references, e.g. [B]. 

In this note we consider the following starting configuration: the cell (0, 0) 
has n particles for some n > 4, every other cell has 2 particles (we call this 
ground level 2). It was shown by Fey-den Boer and Redig j4j in 2008 that this 
model terminates for any n. Define the toppled cluster to be the set of cells 
that were toppled at least once during the stabilisation. Similarly, the visited 
cluster are the cells that were either toppled or received a particle during the 
stabilisation. The toppled and the visited clusters are closely related — it is 
easy to see that the latter is obtained by taking the former together with its 
outer boundary. 

An important consequence of the abelian property is the so-called mono- 
tonicity of the model, which means, adding more particles to the initial 
configuration can only increase the toppled, and therefore also the visited, 
cluster. This naturally looking statement needs to be proved rigorously, and 
this is done by choosing an appropriate order of topplings. 'Ignore' the added 
particles and run the stabilisation as if they were non-existent. We obtain 
the 'old' toppled cluster. Now add the new particles and stabilise the ob- 
tained configuration. The toppled cluster can only increase; that proves the 
monotonicity. 

Define the radius of a set of cells to be sup(^ y)^,^ max(|a;|, |?/|) -|- 1 if 
is non-empty and if is empty. We are interested in the shape and the 
radius of the toppled cluster for our family of starting configurations. The 
shape is determined by the following theorem, proved by Fey-den Boer and 
Redig jl]. 

Theorem 1. For every positive integer n > 4, the starting configuration 
with n particles at the origin and ground level 2 terminates, and the toppled 
cluster is a square centred at the origin. 

The result we prove here concerns the radius of Tn- The best previously 
known bounds of ^/n < r < n/4 for the radius of the cluster were essen- 
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tially trivial. It was conjectured by Fey-den Boer and Redig |1] in 2008 and 
suggested by computer simulations that ^Jn is the right order of magnitude. 
This was proved by by Fey, Levine and Peres [3] in 2010. Our aim here is to 
give a much simpler proof. 

Theorem 2. The toppled cluster of the ahelian sandpile for n particles at 
the origin and ground level 2 is a square, centred at the origin, of radius at 
most 4(1 + o{l))^Jn. 

The main difficulty in proving Theorem [2] lies in choosing the right ap- 
proach. Trying to calculate the radius, let alone the number of particles in 
any given cell exactly, as an explicit function of n, seems to be a very difficult, 
if not a hopeless task. Even though the final configuration exhibits a certain 
pattern that recurs for most values of n (see 0,11]), at the moment of writing 
we are very far away from understanding it. Therefore, in order to prove 
Theorem [2] we had to find a method of estimating the radius asymptotically, 
without referring to the exact numbers of particles in the final configuration. 
The idea of combining the statements of Lemma [1] and Lemma [2] is harder 
to discover than it might seem at the first glance. 

Proof. We choose the following order of topplings: first ignore 2 particles 
everywhere and run the stabilisation with n — 2 particles at the origin on 
ground level 0. Lemma [T] shows that the toppled cluster of this stabilisation 
has a radius of order ^/n. After that Lemma [2] makes sure that adding one 
particle everywhere expands the toppled cluster by at most a factor of 2. 
Adding a particle everywhere twice, therefore, can increase the cluster by at 
most a factor of 4, which gives the desired bound. 

Lemma 1. The radius of the toppled cluster of the ahelian sandpile for n 
particles at the origin and ground level is at most (1 -|- o{l))^Jn. 

Proof. Consider two neighbouring cells of the toppled cluster. Let us call 
them a and h. Without loss of generality we may assume, that the last time a 
or h was toppled, it was h. Then, by the toppling rule, a must contain at least 
1 particle in the final configuration. By a result of Fey-den Boer and Redig [1] , 
the proof of which we do not present here, the toppled cluster contains a 
diamond of radius r, i.e. Dy = {{x,y) : \x\ + \y\ < r — 1}. Subdivide Dr 
into 1 X 2-rectangles, each consisting of two cells; cells that may remain 
are negligible. We obtain (1 -|- o(l))r^ such 1 x 2-rectangles. In the final 
configuration each of them must contain at least one particle. Since the 
ground level was 0, all these particles must come from the origin. This shows 
that n>{l + o(l))r2, and thus r < {1 + o{l))^/n. □ 
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Lemma 2. The square Sr filled with 4 particles on ground level 2 grows by 
at most a factor of 2, i.e. the toppled cluster of its stabilisation has a radius 
of at most 2r. 

Proof. More generally, we show that if Sj.^ is filled with 4's and \ Sj.^ is 
filled with 3's on ground level 2, then the toppled cluster of the stabilisation 
is contained in S^i+ra- We do it by induction on ri for a fixed ri + r2. If 
Ti = 0, then nothing is to prove. Now suppose, ri = k and we know that the 
statement is true for smaller values of vi. Topple every cell of once; this 
is possible by choosing an appropriate order of topplings, as in the proof of 
Theorem [TJ Analogously, the resulting configuration consists of full of 4's 
inside Sr2-i full of 3's, surrounded by a frame of 2's and I's and further 3's 
outside. Repeat the procedure of toppling the square filled with 3's and 4's 
r2 — k times, until we obtain Sk full of 4's, surrounded by a frame of I's and 
2's and stable cells outside. Next, topple every cell in 5*^ once. Note that 
so far the toppled cluster is contained in S'^^+ra- The obtained configuration 
contains 5*^-1 full of 4's, all other cells being stable, and no cell outside Sr2+i 
has 3 particles. Therefore, this configuration has fewer particles in each cell 
than the square of 4's of radius k — 1 inside the square of 3's of radius r2 + 1, 
whose toppled cluster is, by induction hypothesis, contained in S'^.^+ra- Hence, 
by monotonicity, the same holds for the original configuration. □ 

Now we are ready to prove Theorem [2l We start out with n — 2 particles 
at the origin and everywhere else, ignoring 2 particles in each cell; run the 
sandpile on this configuration. By Lemma [1] we obtain a toppled cluster of 
radius at most ri = {l + o{l))^/n. Now we add one particle everywhere. The 
obtained configuration has fewer particles in each cell than the square of 4's 
of size ri on ground level 2, therefore, by Lemma [2] and monotonicity the 
toppled cluster grows by a factor of at most 2. Now add the second particle 
everywhere. The obtained configuration has at most as many particles in 
each cell as the square of 4's of radius 2ri on ground level 2. Hence, by 
monotonicity, the toppled cluster grows again by a factor of at most 2. The 
radius of the resulting final configuration is thereby at most 4ri = 4(1 + 
o(l))v^- 

□ 

Note that Theorem [2] and its proof can be transferred to d dimensions 
(for ground level 2d — 2), showing that in Z'' the radius of the toppled cluster 
is at most CdU^^^, where q is a constant depending solely on d. 
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